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Abstract. We present the first measurement of temperature and polarization angular power spectra of the diffuse 
emission of Galactic dust at 353 GHz as seen by Archeops on 20% of the sky. The temperature angular power 
spectrum is compatible with that provided by the extrapolation to 353 GHz of IRAS and DIRBE maps using 
IFinkbeiner et al. 19991 model number 8. For Galactic latitudes |6| > 5 deg we report a 4 sigma detection of large 
scale (3 < I < 8) temperature-polarization cross- correlation (£ + l)Cf E /2n = 76 ±21 fiK R j 2 and set upper limits 
to the E and B modes at 11 //Krj 2 . For Galactic latitudes |6| > 10 deg, on the same angular scales, we report a 2 
sigma detection of temperature-polarization cross-correlation (£ + l)Cf E /2n = 24 ± 13 /iKrj 2 . These results are 
then extrapolated to 100 GHz to estimate the contamination in CMB measurements by polarized diffuse Galactic 
dust emission. The TE signal is then 1.7 ± 0.5 and 0.5 ± 0.3 /iKq MB for |6| > 5 and 10 deg. respectively. The 
upper limit on E and B becomes 0.2 /jKq MB (2a). If polarized dust emission at higher Galactic latitude cuts is 
similar to the one we report here, then dust polarized radiation will be a major foreground for determining the 
polarization power spectra of the CMB at high frequencies above 100 GHz. 

Key words. Cosmic Microwave Background - Cosmology - Observations - Submillimetre - Polarization - Dust - 
Extinction - Foreground 



1. Introduction 

The Cosmic Microwave Background (CMB) is now con- 
sidered as one of the most sensitive probes to the physics 
of the early Universe. A great number of experiments 
have measured its temperature anisotropy power spec- 
trum over a wide range of angular scales (for a review, 
see [Wang et al. 2003| l until WMAP recently gave a cos- 
mic variance limited estimate up to the rise of the second 
acoustic peak l|Bennett et al. 20 03a). 

CMB polarization provides a wealth of complemen- 
tary information. First, it allows to break degenera- 
cies that remain in the determination of cosmologi- 
cal parameters with temperature anisotropy data alone 
(see e.g. |Zaldarriaga et al. 1997| ). Perhaps more impor- 
tantly, it sheds light directly on inflation through the B 
mode that is generated only by primordial gravity waves 
produced during that era ( |Seljak fc Zaldarriaga 19971 
IKamionkowski et al. 1997)) . Distortions of the E mode by 
large scales structures weak lensing also induce a Bwl 
signal, distinguishable from the primordial B by its non- 
gaussianity. This Bwl component in turn provides use- 
ful information about the dark matter distribution (e.g. 
IBenabed et al. 20T)T)l and the mass of neutrinos at the pre- 
cision of 0.04 cV for an experiment e.g. 20 times more 
sensitive than Planck fKaplingh at et al. 2003| ). 

CMB polarization is however 2 to 5 orders of magni- 
tude smaller than temperature anisotropies and therefore 
still remains to be accurately measured. It is now becom- 
ing accessible thanks to improved instrumental sensitivi- 
ties. The first detection of the E mode has been reported 
by DASI IjKovac et al. 2002|> , This result has been con- 
firmed by the same team l|Leitch et al. 200411 . CAP MAP 
HBarkats et al. 2004)1 and CBI IjReadhead et al. 2004)1 . 
WMAP has also provided an estimate of the TE corre- 
lation spectrum jKogut et al. 20 03 ) fully compatible with 
an inflationary scenario. 

Instrumental sensitivity is not the only issue in the de- 
termination of the power spectrum of the polarized CMB 
anisotropies. Other astrophysical components also con- 
tribute to the sky brightness and polarization at the wave- 
lengths of interest and must be subtracted. These fore- 
grounds, mainly dominated by the diffuse Galactic emis- 
sion and the one of extragalactic point sources, are often 
not well constrained or even not experimentally measured 
in the case of polarization. 

Ground based experiments such as those mentionned 
above observe small regions of the sky. They can then 
choose them where foregrounds are weak and are less 
prone to contamination by Galactic emission. It is not 
the case for full sky experiments such as WMAP and 
Planck. At the Planck-HFI frequencies and for future 
bolometer experiments, the dominant component is the 
radiation from Galactic Interstellar Dust (ISD). The sub- 
Send offprint requests to: reprints@archeops.org 
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millimetre and millimetre (hereafter submm) intensity of 
the ISD emission can be inferred from IRAS and COBE- 
DIRBE extrapolations (e.g. Finkbcinc r et al. 19991 here- 
after FDS) and has been measured on large scales by 
COBE-FIRAS HReach et al. 19951 |Boulanger et al. 19961 
|Lagache et al. 1998| ). However, little is known about ISD 
polarization emission on scales larger than 10 arcmin, 

1. e. that are the most relevant for current CMB studies. 
Indeed, ground-based observations of submm ISD polar- 
ization are concentrated on high angular resolution (ar- 
cminute scale) star formation regions. Indirect evidence 
for large scale polarization comes from the polarization of 
starlight in extinction (ISerkowski et al. 1975)1 . but as re- 
viewed by Goodman (1996), these measurements of back- 
ground starlight polarization lead to ambiguous interpre- 
tation. In particular, the visible data are biased by low col- 
umn density lines of sight and do not fairly sample more 
heavily reddened ones. Direct submm measurements are 
therefore highly required both for Galactic studies of the 
large scale coherence of the magnetic field and in the field 
of CMB polarization, but are rather challenging as they 
require sensitivities comparable to those of CMB studies. 

Recently, IBenoit et al. 200"4al have reported the first 
measurement of the submillimetre diffuse polarized emis- 
sion by interstellar dust in the vicinity of the Galactic 
plane using the Archeops experiment. They show that 
the Galactic plane is significantly polarized at the 3-5% 
level and that dense clouds can be polarized up to 10% 
or more. This indicates that the dust intrinsic radiation is 
highly polarized and that the grain alignment mechanism 
is very effective. Considered with the possible large scale 
coherence of the polarization orientation, it shows that the 
dust polarized emission could be an important foreground 
for CMB polarization studies, especially on large angular- 
scales. 

Here, we wish to give a first answer to this question 
with the evaluation of the dust polarization power spec- 
tra away from the Galactic plane and in diffuse regions, 
on angular scales ranging from I = 3 to £ = 70. Section [3 
briefly introduces Archeops and its polarization capabili- 
ties. Section |3 presents the processing applied to the data 
and Sect. 2] describes the evaluation of the polarized angu- 
lar power spectra. Section^lpresents our main results. We 
conclude in Sect. with the extrapolation of our results 
to lower frequencies where the CMB dominates, to give an 
estimate of the dust contamination in the measurements 
of the CMB polarization power spectra. 

2. The Archeops instrument 

Archeops 1 is a balloon borne bolometer experiment that 
aimed at measuring the CMB temperature anisotropy 
over large and small angular scales. It provided the 
first determination of the Ci spectrum from the COBE 
multipoles IjSmoot et al. 19 92) to the first acoustic peak 
l)Benoit et al. 2003a)) from which it gave a precise deter- 

1 http://www.archeops.org 
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urination of the main cosmological parameters, such as 
the total density of the Universe and the baryon frac- 
tion (JBenoit et al. 2003b|) . Archeops was also designed as 
a test bed for Planck-HFI and therefore shared the same 
technological design: a Gregorian off-axis telescope with 
a 1.5 m primary mirror, bolometers operating at common 
frequencies (143, 217, 353 and 545 GHz) cooled down at 
100 mK by a 3 He/ 4 He dilution designed to work at zero 
gravity and similar scanning strategy. A detailed descrip- 
tion of the instrument and its performances can be found 
in IjBenoit et al. 2002|> and IjBenoit et al. 2003b|> . 

At 353 GHz where dust thermal radiation is dominant, 
Archeops has 6 bolometers mounted in 3 Ortho Mode 
Transducer 2 (hereafter OMT) pairs that are sensitive to 
polarization in order to study the properties of the polar- 
ized dust diffuse radiation. The three OMTs are oriented 
at 60 degrees from each other to enable the full recovery of 
the Q and U Stokes parameters and to minimize the corre- 
lations in their determination. Archeops was launched on 
February 7th, 2002, from the CNES/Swedish facility of 
Esrange, near Kiruna (Sweden). The flight brought about 
12 hours of high quality night data. 

3. Data processing and map making 

A detailed description of the data processing and the po- 
larization map making is given in (Bcn oit et al. 2003b) 
and l|Benoit et al. 2004a|> . 

The main steps on the Archeops data processing are 
the followings. First, the reconstruction of the point- 
ing during flight, with rms error better than 1 arcmin 
HBourrachot 2004jl . is performed using the data from a 
bore-sight mounted optical star sensor aligned to each 
photometer using Jupiter observations. Second, corrupted 
data (including glitches and bursts of noise) in the Time 
Ordered Information (TOI), representing less than 1.5%, 
are flagged and not considered in the following process- 
ing. Third, low frequency drifts on the data, generally 
correlated to house-keeping information (altitude, atti- 
tude, cryostat temperatures, the CMB dipole) are re- 
moved using the latter as templates. Fourth, high fre- 
quency decorrelation is performed in few chosen frequency 
intervals of ~ 1 Hz width to remove non-stationary high- 
frequency noise. Sixth, the corrected timelines are then 
deconvolved from the bolometer time constant and the 
flagged corrupted data are replaced by a constrained real- 
ization of noise. Finally, low frequency atmospheric resid- 
uals and noise are subtracted using a destriping proce- 
dure which preserves the sky signal to better than 2% on 
large angular scales (Bou rrachot 2 004). To improve the 
quality of the removal of atmospheric residuals we have 
also performed a component separation analysis in the 
time ordered data using the SMICA-MCMD algorithm 
IjDelabrouille et al. 2003|l over all the Archeops channels. 
From this analysis we have constructed a template of the 



atmospheric contribution to the Archeops data which has 
been fit and removed from each of 353 GHz bolometers 
preserving the dust emission to better than 5%. 

The six bolometers are cross-calibrated as discussed 
on Sect. 14.3.21 The absolute calibration is obtained from 
a correlation between the Galactic latitude profiles from 
FIRAS "dust spectrum" maps and those of Archeops. It 
has an absolute accuracy of about 12%. This affects only 
the absolute values of /, Q, U but neither the degree of 
polarization nor its orientation. A detailed description of 
the calibration is given in IBenoit et al. 2003bl 

The data processed and calibrated as above have 
been projected into polarization maps using the algo- 
rithm described in [Ponthieu 2003 The maps produced 
are shown on Figs. Q [3 El A significant polarization 
level in the Galactic plane was first reported by Archeops 
( Bcn oit et al. 2004a)) for regions with longitudes between 
100 and 120 degrees and between 180 and 200 degrees. The 
new data processing allows us to reconstruct polarization 
on all the Galactic plane observed by Archeops and to fill 
the gap between 120 and 180 degrees. 

The noise correlation matrix on those maps have been 
computed using pure noise simulations which are de- 
scribed in Sect. 14.21 



I 
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Fig. 1. Total unpolarized intensity / measured by 
Archeops at 353 GHz. Map centered on Galactic longi- 
tude I = 120 in Galactic coordinates. The pixel size is 
27 arcmin smoothed by a 2 deg beam FWHM to match 
the 1.88 deg pixel size (HEALPix parameter n S id e = 32) 
used throughout the analysis. Grid lines are spaced by 
20 deg. 



2 Planck-HFI has since changed for the Polarized Sensitive 
Bolometer technology to measure polarization of radiation 
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Fig. 2. Stokes parameter Q measured by Archeops at 
353 GHz. Map centered on Galactic longitude I = 120 in 
Galactic coordinates. The pixel size is 27 arcmin smoothed 
by a 2 deg beam FWHM. Grid lines are spaced by 20 deg. 
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Fig. 3. Stokes parameter U measured by Archeops at 
353 GHz. Map centered on Galactic longitude I = 120 in 
Galactic coordinates. The pixel size is 27 arcmin smoothed 
by a 2 deg beam FHWM. Grid lines are spaced by 20 deg. 



4. Evaluation of the polarization power spectra 

4.1. Formalism 

For a direction of observation n, we define the Stokes pa- 
rameters /, Q and U in the tangential plane with respect 
to (— e^e^,). The angle of the polarization is oriented 
from the North Galactic pole through East to the South 
Galactic pole (counterclockwise). 



For statistical analysis, the use of E and B 
QSeljak fc Zaldarriaga 1997| ) is however more suitable be- 
cause these quantities are scalar and independent of the 
coordinate system. We here estimate these quantities and 
the correlation between E and T using the method de- 
scribed in IChon et al. 20041 The spin-2 nature of Stokes 
parameters leads us to define 



P = Q + iU. 



(1) 



When the polarized two point correlation functions are 
estimated in real space between two directions h\ and h 2 , 
one has to rotate P by angles a\ and a 2 respectively to 
align the axes defining Q for each direction of observation 
with the geodesic connecting hi and h 2 such that 



P{h 2 ) 



„2ia 



'P(hi) 
e 2ia 'P(h 2 ). 



(2) 
(3) 



We then define the correlation functions 



= (P(hi)P(n 2 )) 
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(6) 



where hi -hi — cos 9 and d! mm , are the reduced Wigner 
Z)-matrices. The Cgs angular power spectra are then ob- 
tained by the following integration: 



r<EE r<B o- n EB 

cf E + cf 
cJ E + icT B 



2vr / L (d)4_ 2 (0)d cos 0, (7) 



2ir J J + (8)4 2 (6)d cos 0, (8) 
2tt / £ x (0)d e 2Q d cos 6 . (9) 



We estimate those quantities with the software 
SpicePol l|Ghon et al. 2004jl that uses the HEALPix pack- 
age (Gorski et al. 1998jl to compute the pseudo-C.es from 
the raw maps : Cf E , Cf E , Cf B , from which we obtain 
an estimate of the signal plus noise angular power spec- 
tra. The noise contribution is estimated through Monte- 
Carlo simulations (Sect. I4.2|l and subtracted, producing 
estimates of the angular polarized power spectra. Because 
of approximations used by SpicePol it can not produce 
estimates of Cj B and Cf B . 

The total uncertainty on the determination of the 
power spectrum Cj T includes both the variance of the 
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noise and the variance of the signal itself. In the case of 
Galactic dust, there is no cosmic variance stricto sensu. 
Indeed, the dust distribution in the sky and its proper- 
ties are not a particular realization of a Gaussian random 
field like in the case of the CMB anisotropies. Therefore, 
we compute the variance of the polarized power spectra 
considering only the noise-noise and the signal-noise con- 
tributions, and not the signal-signal covariance (sample 
variance) as follows 



2/ n XX\ 



J2(riTE\ 



(2/ + l)f sky 
1 

{21 + l)f sk y 



N, 



XX 2 



(10) 



\ N JE2 + N TT C EE 



+N EE Cj T + NJ t N ee \ (11) 



where X stands for T, E or B, N is the noise and 
fsky is the fraction of the sky taken for the analysis. The 
dominant term in Eq. Ijllfl is the product of the covariance 
of temperature by the polarization noise N EE Cj T . These 
relations are however only handy approximations to the 
true uncertainties and need an empirical adjustment of the 
fsky parameter that is performed using the simulations 
described in Sect. 14.21 These simulations are also used to 
compute the noise power spectra. 

The temperature power spectrum Cj T and its 
error bars are computed using the Xspect method 
IjTristram et al. 2*0 04) which uses 15 cross power spectra 
from the six detectors and no auto power spectrum in or- 
der to avoid corrections induced by the detector noise. 

4.2. Instrumental noise 

The noise power spectrum of the bolometers is es- 
timated from a four step process already applied in 
Benoit et al (2003a) in the context of the Cj T evaluation. 
First, the Galactic plane region is masked and interpo- 
lated in the time domain with slowly varying functions. 
Secondly, these timelines are projected onto maps that 
are deprojected to obtain a second timeline with a higher 
signal to noise ratio. We then subtract this second time- 
line from the original one to obtain a noise dominated 
timeline and compute its time domain power spectrum. 
Simulations including realistic noise and Galactic contam- 
ination show that this process allows the recovery of the 
true noise power spectrum at the 5% level. 

From these time domain noise power spectra, we gen- 
erate noise timelines for the six polarized bolometers and 
project them onto maps in the same way as the real data. 
The same statistical analysis as the one applied to the data 
is performed on 250 of such noise maps to have a good es- 
timate of the noise angular power spectra. These power 
spectra are then subtracted from the Archeops polarized 
angular power spectra in order to correct them from the 
noise bias. Finally, we compute the uncertainties on the 
polarization power spectra using equations flH| and lfTT|) . 



4.3. Systematic effects 

Three main sources of systematic effects are likely to af- 
fect the evaluation of the polarization correlations: the 
filtering applied to the data, the uncertainty on the cross- 
calibration between the detectors and the uncertainty on 
the knowledge of the exact orientation of the polarizers on 
the focal plane. We address these three issues separately 



4.3.1. Filtering 

A bandpass filtering is applied to the data. The lower fre- 
quency bound is 0.033 Hz, which corresponds to the first 
harmonic f sp i n of the rotation of the gondola (2 rpm). 
Because of the spinning of the instrument, few physical 
components at frequencies below f sp i n in the timeline are 
projected on the map. The high pass therefore ensures 
that no physically irrelevant and dipole-like components 
remain in the timelines. A low pass filtering is then applied 
to the timelines at 38 Hz to remove high frequency noise. 
To correct from this effect on the angular power spectra 
we have computed an effective filtering function Fi in the 
multipole space from simulations of temperature Gaussian 
fields with a flat power spectrum. The maximum correc- 
tion is at low £ and is less than 2%. 

We have also checked that the filtering did not induce 
any spurious polarization such as leaks from total intensity 
into polarization. For this we have applied the bandpass 
filter to simulated timelines, deprojected from the FDS 
template at 353GHz, for the six bolometers involved in 
the determination of the angular power spectra. Then, we 
have reconstructed /, Q and U maps for which we have 
extracted the temperature and polarization power spec- 
tra. No spurious polarization was produced at the level of 
0.1%. 



4.3.2. Cross calibration 

Stokes parameters are mainly estimated from the differ- 
ences of the outputs of bolometers that measure orthog- 
onal polarization states. An error in the cross calibration 
between detectors generates a systematic leak of total in- 
tensity into polarization. The cross calibration of Archeops 
channels is described in details in lBenoit et al. 2004al and 
relies on the scaling of Galactic intensity profiles computed 
from each bolometer. The cross calibration coefficients are 
then determined at the 2% level. In order to give a con- 
servative upper limit to the effect of this uncertainty on 
the angular power spectra, we forced the cross calibration 
coefficients to values which deviate by 2% and maximized 
the relative difference between two orthogonal photome- 
ters. The angular power spectra were then estimated with 
this new set of coefficients and the error bars were derived 
from the standard deviation of 200 simulations, assuming 
a symmetric uncertainty. This uncertainty is about 5% of 
the statistical error bars. 
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4.3.3. Polarizer relative orientation 

The accurate positioning of the polarizers in the fo- 
cal plane is made difficult by the complexity of the in- 
strument. We therefore checked that they were correctly 
placed during the pre-fhght ground calibrations and found 
that they were indeed in their nominal configuration with 
a la uncertainty of 3 degrees. A mismatch between the 
assumed and the real orientations of the polarizers gener- 
ates a relative error of a few percents on the reconstructed 
Stokes parameters (Kapla n et al. 2001] ). In order to esti- 
mate the error induced by the uncertainty of the knowl- 
edge of the accurate positionning of the polarizers, we have 
performed simulations in which these angles are forced to 
random values different from those used for the /, Q, U 
reconstruction. The resulting uncertainty is of the same 
order as that of the cross calibration uncertainty, that is 
to say about 5% of the statistical error bars. 

5. Results 

We estimate the angular power spectra as a function of 
cuts of the data in Galactic latitude in order to remove the 
effects of regions with the strongest dust emission along 
the Galactic plane. Figure 0] gives the power spectra for 
\b\ > 5, \b\ > 10, and \b\ > 20. The spectrum for Cj T oscil- 
lates as a function of £, with the amplitude of oscillation 
increasing with decreasing £. These features are consis- 
tent with a Galactic origin for the signal as discussed in 
Sect. EXT] 

The TE Spectrum: on large angular scales 3 < £ < 8 
there is a ~ 4cr detection for \b\ > 5 with a magntidue of 
(£ + l)Cf E /2ir = 76 ± 21 ^Krj 2 and a ~ 2a detection 
for |6| > 10 with a magnitude of (£ + l)Cj E /2tt = 24 ± 
13 /xKrj . For |6| > 5 there is also a ~ 2 a detection for 
18 < £ < 23, corresponding to a peak in the temperature 
spectrum at the same I bin. 

The EE Spectrum: on large angular scales 3 < £ < 8 
there is a ~ ha detection for |6| > 5 with a magnitude 
of (£ + l)Cf E /2n = 7.5 ± 1.5 /iK RJ 2 and a - 2 a for 
18 < £ < 23. Otherwise, the power in the EE spectra is 
consistent with zero for all £'s for |6| > 10 and for |fe| > 20 
and at 20 < £ < 70 for all latitude cuts. Similar to the 
TE spectrum, the power measured for EE on all scales 
decreases with increasing latitude, which is consistent with 
a Galactic origin for the signal. 

The BB Spectrum: there are detections of power in 
the first few £ bins and the spectrum is consistent with 
no power for £ > 19. The power at the low £ bins does 
not depend on Galactic latitude cut. We argue below that 
models of polarization from dust emission predict a de- 
crease of power in both the EE and BB spectra with 
Galactic latitude cut and therefore the BB signal at low 
£ is probably spurious. Unlike the CMB for which E and 
B have different physical sources, dust is expected to pro- 
duce comparable amounts of E and B. Assuming a similar 
level of E and B and interpreting the data in terms of up- 
per limits we provide a common 2a upper limit to both 



modes at (£ + l)Cf E ' BB /2tt < 11 ^Krt 2 at 353 GHz for 
|6| > 5. 

Since the Archeops data probe the dust diffuse emis- 
sion and not local clouds in the Galactic plane it is likely 
that the orientation and coherence of the magnetic field is 
similar for latitudes larger than 5 and larger than 10 deg. 
Therefore, the more dust there is along the line of sight, 
the more intense the emission should be, and the more 
power we expect in the temperature and polarization 
spectra. According to these arguments there should be 
a monotonous decrease of power in the spectra with in- 
creasing latitude cuts, as observed for the TT, TE and 
EE data. 

5.1. Comparison with models 

In this section we compare our results to models of diffuse 
Galactic dust emission. First, we compare the TT power 
spectrum to those expected using a Galactic cosecant-law 
model and using the FDS model. Second, we obtain an al- 
ternative estimate of the TE spectrum by cross-correlating 
the FDS template with the Archeops Q and U maps. 
Finally, we compare the measured spectra to those cal- 
culated on the basis of a simplified physical model of po- 
larized emission from dust. 

5.1.1. Pure cosecant-law Galactic emission 

Due to the disk geometry of the Galaxy, the integrated 
emission along a given line of sight increases as the abso- 
lute value of the latitude decreases. This is well approx- 
imated by a cosecant law of the form 1(b) cx l/sin(6). 
When this Galactic contribution is subtracted, the angu- 
lar power spectrum of the remaining diffuse dust has the 
form Ct cx £~ 3 IjGautier et al. 19921 |Wright 1998| ). 

We choose to leave the Galactic contribution in our 
data such that we can assess its magnitude and its poten- 
tial contamination of CMB polarization experiments, (see 
Sect. [BJ. In order to estimate the influence of a cosecant 
law component in our data, we simulate such an emis- 
sion analytically and with an amplitude compatible with 
the FDS template. The results are presented on Fig. [5] 
where we plot the TT angular power spectrum at differ- 
ent Galactic latitude cuts for the Archeops data, for the 
cosecant-law model and for the FDS template. There is a 
qualitative agreement between the oscillation pattern of 
the Archeops data for low Galactic latitude cuts and the 
pattern present for the cosecant-law model and for the 
FDS template. The agreement suggests that the cosecant- 
law emission dominates the observed Galactic dust emis- 
sion at large angular scales. 

5.1.2. FDS-Archeops correlation 

It is interesting to correlate the Archeops Q and U 
data with available dust intensity maps. We quantify 
the correlation by calculating the cross-spectrum be- 
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Fig. 4. Clockwise from top left: power spectra Cj T , Cf E 7 Cf B and Cf E computed from the 353 GHz Archeops data 
for three different Galactic cuts \b\ > 5, 10, and 20 deg. At low I the power in Cf T , Cf E , and Cf E decreases with 
increasing |6|, as would be expected from a Galactic signal. Since the power in Cf B does not change with |6| its source 
is probably not astrophysical. ASCII files of these data can be obtained at http://www.archeops.org/info_polar.html 



tween the two data sets. This idea has already been 
succesfully applied to temperature anisotropy data sets 
IjAbroe et al 20041 ITristram et al. 20041 For dust inten- 
sity we use the FDS template, based on IRAS 100 /im 
maps and extrapolated to 353 GHz ("model number 8") 
HFinkbeiner et al. 1999J1 . We found that this model fits 
the Archeops temperature data at 353 GHz with frac- 
tional deviations of less than 20%. Figure |H1 shows that 
the "FDS x Archeops" cross-spectrum (in blue) and the 
"Archeops x Archeops" TE power spectrum (in red) are 
consistent within the error bars. This result provides ad- 
ditional confidence that the detection of power in the 
Archeops x Archeops TE spectrum at low I is due to po- 
larized dust and is not spurious. We note that the cross- 
spectrum with the FDS map avoids correlated noise that 
is inherent in the Archeops x Archeops spectrum and thus 
the agreement between the two results suggests that cor- 
related noise did not induce spurious results. 



5.1.3. Comparison to a simple physical model of 
Galactic dust polarization 

In a given direction of observation, the measured Stokes 
parameters are the result of the integral along the line 
of sight of the local Stokes parameters. These, in turn, 
depend on the local alignment of dust grains with the 
magnetic field and on the intrinsic degree of polariza- 
tion. Precise modelling of the diffuse emission due to 
dust and its polarized angular power spectra is a com- 
plex problem. Previous work has shown that on large 
scales, the alignement of dust grains was compatible 
with a Galactic magnetic field aligned with the spiral 
arms (|Fosalba et al. 20021 IBenoit et al. 2004al and refer- 
ences therein). However, such work relate to data at low 
Galactic latitudes \b\ < 10 deg, and the extrapolation at 
higher latitudes yet remains uncertain due to lack of data. 
We here use a toy model of extrapolation at high lati- 
tudes of these polarization properties that gives qualita- 
tive agreement with the Archeops data. 



Since we are here looking at latitudes away from the 
Galactic plane, most of the dust that is being probed is in 
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Fig. 5. Temperature angular power spectrum of Archeops data (upper left diagram), of a pure Galactic cosecant 
emission law simulation (upper right diagram) and of a FDS template (bottom diagram) computed using Xspect for 
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Fig. 6. TE angular power spectra for \b\ > 5 either with Archeops data alone (red, circles) or with Archeops Q and U 
data and / from the FDS template (blue, diamonds) on the Archeops region of observation extrapolated to 353 GHz 
l|Finkbeiner et al. 1999j> . 



our vicinity. Typically, if we take the Galactic disk to be 
200 pc thick, a line of sight at b = 5 deg exits the plane 
at ~ 1.1 kpc from the observer. This is small compared 
to the size of the spiral arm in which we are located, and 



we therefore assume that the large scale magnetic field 
component aligned with the local spiral arm is along a 
constant direction on this portion of space. A realistic 
model would of course consider turbulence, but we are 
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only interested here in a first order approximation. Based 
on lBenoit et al. 2004al who showed that dust diffuse emis- 
sion in the vicinity of the Galactic plane was polarized at 
the 5% level, and that some dense clouds were polarized 
at more than 10%, we present in Fig. [7| results with three 
assumptions for the level of polarization of Galactic dust: 
p(x) — 5, 10, and 15%. We have only considered latitude 
cuts of \b\ > 5 deg. Given the simplicity of the model, the 
relative agreement between the data and the shape and 
amplitude of the models is encouraging. 

6. Dust contamination in CMB polarization 
estimates 

Our results can be used to estimate the contamination of 
dust to CMB polarized angular power spectra. For this we 
need to extrapolate our measurements at 353 GHz to lower 
frequencies where the CMB is usually measured since it is 
where it is more intense. 

Note that this extrapolation concerns large dust grains 
only. It does not take into account the so-called "anoma- 
lous emission", since (1) the latter does not contribute sig- 
nificantly at 353 GHz and (2) it is likely not to be polarized 
(e.g. ILazarian fc Prunet 2001(1 . Dust thermal radiation is 
the dominant foreground at frequencies above ~ 90 GHz. 
At high frequencies it has long been described by a single 
grey body at a temperature of 17.5 K with a v 2 emissiv- 
ity (Boulanger et al. 1996). At low frequencies (~ 300 to 
~ 100 GHz), the spectrum is not well constrained due to 
the lack of accurate measurements. We need to use the 
FIRAS and WMAP data jointly. The FIRAS data are 
used to compute the intensity value at 353 GHz and the 
WMAP data give the intensity value at 94 GHz. This then 
allows to derive the dust radiation spectral index. 

To determine these dust intensities we use the method 
described in Lagache (2003) and derive the spectrum of 
the dust whose emission is correlated with the HI gas, 
using FIRAS and WMAP data in the Archeops region. 
We obtain a ratio J„ (353)//,, (94) = 134 corresponding 
to a spectral index of 1.7 that we use in the following. 
Note that this value is very stable when different (large) 
regions of the sky are averaged. This is also the same value 
as the one obtained by Finkbeiner et al (1999) in this 
frequency range (see discussion in lFinkbeiner 2004|l . Using 
this spectral index of 1.7, we find for the largest angular 
scales (3 < £ < 8), {£ + l)Cj E /2tt = 1.7 ± 0.5 and 0.5 ± 
0.3 fj,K^ MB for \b\ > 5 and 10 deg. respectively at 100 GHz. 

These results suggest that dust may be a very signif- 
icant foreground for measurements of the CMB polariza- 
tion angular power spectra, particularly if they include 
Galactic latitudes below |6| < 10 degrees. However, these 
measurements must be interpreted with care when con- 
sidering the entire sky. Our data includes only 20% of the 
sky and generalization to the entire sky is questionable 
because of the complexity of the statistical properties of 
dust. 

It is interesting to consider the relation between our 
results and those of the WMAP team. By combining data 



from five frequency bands at lower frequencies (23, 33, 
41, 61 and 94 GHz) and for the entire sky the WMAP 
team has computed a value for the TE spectrum at 
3 < t < 8 of 1.72 ± 0.50 Ai K cMB- Accurate extrapola- 
tion of the Archeops results into the WMAP data is com- 
plicated and is beyond the scope of this paper. (For ex- 
ample, it would require a detailed knowledge of the com- 
ponent separation methods, the relative weighting of the 
frequency bands and the shape of the beam). However, it 
is easy to illustrate that combining data at lower frequen- 
cies reduces the effects of dust substantially. We extrap- 
olate our results to the WMAP frequency bands with a 
constant emissivity spectral index of 1.7 in antenna tem- 
perature. We use an equal weight per frequency and use 
only the area of the sky that overlaps the Archeops data 
and WMAP's Kp2 mask. The result of this extrapolation 
is presented on Fig. [S] On scales of 3 < £ < 8 we find 
(£ + l)Cf E /2ir = 0.17 ± 0.06 fiK% MB , which is about a 
factor of 10 smaller than both our estimate at 100 GHz, 
and WMAP's reported result. 

When extrapolated to 100 GHz the upper limits that 
we found on the E and B modes for |6| > 5 become 
0.2 /iK^ MB . We note that the CMB B mode is expected to 
be at most ~ 10~ 3 /iK^ MB at I ~ 5, if the reionization op- 
tical depth is as high as r = 0.17 fSpergel et al. 20031 ) an d 
if the tensor to scalar ratio T/S is the highest compatible 
with current CMB temperature anisotropy measurements. 
If the actual level of polarized dust over most of the sky is 
close to the upper limit we found in this work then a sub- 
traction of a large foreground signal will be necessary even 
at 100 GHz in order to detect the primordial gravitational 
waves on large angular scales. 

7. Conclusions 

In this paper we have presented the first measurements 
of the angular power spectra of the Galactic dust polar- 
ized diffuse emission on approximately 20% of the sky 
at 353 GHz by Archeops. On angular scales 3 < I < 8, 
we obtain a 4 a detection of (I + l)Cj E (dust) / 2n = 
76 ± 21 /zKrj 2 for |6| > 5 deg. On the same angu- 
lar scales and for \b\ > 10, we have a 2 a detection 
(£ + l)Cj E {alust)/2ii = 24 ± 13 ^Krj 2 . This decrease in 
power is expected from the cosecant behaviour of the large 
scale component of the total intensity, which is shown to 
dominate the total intensity angular power spectrum. 

On the same sky coverage and for all angular scales 
3 < I < 70, we set upper limits to the E and B polar- 
ization at (£ + l)C EE ' BB (dust) /2tt < 11 ^K RJ 2 . These 
results have been confirmed by using a template of the 
Galactic dust intensity from Finkbeiner et al (1999) in 
place of the Archeops total intensity map to compute the 
TE angular power spectrum. This spectrum agrees with 
the one derived from Archeops data alone on all angular 
scales within 1 a. 

Furthermore, the high degree of polarization seen in 
the Galactic plane by Benoit et al (2004a) together with a 
simple model of the Galactic magnetic field and the align- 
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Fig. 7. Comparison between Archeops polarized angular power spectra and those of a simple model of the Galactic 
magnetic field for all latitudes \b\ > 5. The large scale magnetic field is assumed to be constant and oriented along 
the local spiral arm. Dust is assumed to be aligned orthogonaly to it. The effective degree of alignment of the grains 
together with the dust intrinsic polarized emissivitiy is assumed to lead to 5 (blue), 10 (red) and 15% (green) effective 
degree of polarization. 



ment of dust grains, leads to estimates that are compatible 
with the data. 

To estimate the contribution of Galactic dust to the 
measurement of polarized CMB anisotropics, we have ex- 
trapolated our results to the reference frequency 100 GHz, 
using a spectral index of 1.7, derived from FIRAS and 
WMAP data on the Archeops sky coverage. The TE 
mode becomes (£ + l)Cf E {dust)/2n = 1.7 ± 0.5 and 
0.5±0.3 mKcmb on 3 < ^ < 8 for |6| > 5 and 10 deg respec- 
tively. The upper limit on the E and B modes becomes 
0.2 fiK^ MB . These values show that even at 100 GHz 
where dust radiation is expected to be lower than the 
CMB, its polarization may be very significant compared 
to the CMB and should be subtracted with care from the 
observations. 

The effects of the polarization of dust are less severe at 
low frequencies. When extrapolating our measurements at 
353 GHz to 20 GHz with the constant spectral index 1.7 
and weighting the 5 WMAP frequencies equally we find 
a level of (I + l)Cj E '{dust) /2tt = 0.17 ± 0.06 A*K^ MB for 
3 < I < 8 on the intersection between WMAP's Kp2 mask 



and Archeops sky coverage. This is about a factor of 10 
smaller then the WMAP team results. 

The high level of polarization measured at 353 GHz 
by Archeops and anticipated at the reference frequency of 
100 GHz, together with the uncertainties on dust spectral 
index and the extrapolation of its statistical properties to 
the whole sky call for further precise studies in order to 
subtract it precisely from CMB data. This is even more 
critical for the detection of the imprint of the primordial 
gravitational waves on the CMB B mode, that is expected 
to be much smaller than the present upper limit on dust 
B polarization. 
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